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Abstract 

We consider linear cosmological perturbations on the background of a D-brane gas in which the 
compact dimensions and the dilaton are stabilized. We focus on long wavelength fluctuations and 
find that there are no instabilities. In particular, the perturbation of the internal space performs 
damped oscillations and decays in time. Therefore, the stabilization mechanism based on D-brane 
gases in string theory remains valid in the presence of linearized inhomogeneities. 
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I. INTRODUCTION 



It is well known that the theory of big-bang nucleosynthesis together with the obser- 
vational abundances of light elements provide a limit to the change in the size of internal 
space following the production of these elements. Therefore, string theory should be able 
to explain why the extra dimensions are stabilized following nucleosynthesis? Intuitively, 
winding branes are expected to resist the expansion and this may offer an explanation. An 
accommodation of this idea in early string theory can be found in the paper of Branden- 
berger and Vafa [1], and recently the stabilization of toroidal extra dimensions by winding 
strings is verified by explicit calculations in [2, 3] (see also [4]). The model of [1] is developed 
in [5] to include other extended objects in string theory and this scenario is now known as 
brane gas cosmology (BGC) (see e.g. [6]- [31] for recent work). 

It was also claimed in [1] that the dimensionality of space-time can be explained by the 
annihilation of winding strings in a three dimensional subspace leading decompactification. 
From the naive intersection probability of hypersurfaces embedded in d-spatial dimensions 
one can infer that d = 2p + 1 is the critical value which allows annihilation of p-branes. 
Remarkably, for strings d = 3 and according to [1] that is why a three dimensional subspace 
is decompactified. As an additional consequence of this argument, in BGC all higher dimen- 
sional branes are expected to be annihilated in the early universe. Although these claims 
are supported numerically in [32], there are also some caveats. It is observed in [15] that the 
number of decompactified dimensions depends on initial conditions and the outcome of three 
large dimensions is argued to be statistically insignificant [25]. In [26] it is shown that the 
interaction rates of strings are not strong enough to sustain thermal equilibrium. Moreover, 
D-brane anti-D-brane annihilation rate turns out to be small compared to the expansion 
rate of a FRW universe including the higher dimensional D-branes [12]. Therefore, it seems 
that the dimensional counting argument of [1] is not sufficient alone for the annihilation 
process and it is plausible to have surviving branes other than strings or membranes. 

In the light of these results, we think that late time models with higher dimensional 
branes should also be studied in the context of BGC. Such models may be useful in solving 
two difficulties with string gases. Firstly, as pointed out in [7], strings may not be able to 
stabilize volume modulus, or further, there may not be any winding strings in the spectrum 
depending on the topology of internal space (consider e.g. a sphere). Secondly, it turns out 
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to be difficult to stabilize the dilaton and keep up the stability of extra dimensions [23], 
which may cause problems in canonical Einstein frame after compactification. 

In [24], we have considered a gas of D6-branes wrapping over the six- dimensional internal 
space and demonstrated dynamical stabilization of the volume modulus and the dilaton. The 
aim of this paper is to consider long wavelength perturbations in this model to see the validity 
of the conclusions under linearized inhomogeneities. For the string gas the perturbation 
analysis has been carried out in [18, 19]. The plan of the paper is as follows. In the next 
section we derive linearized equations for the scalar metric and dilaton perturbations in the 
generalized longitudinal gauge. In section III, we focus on long wavelength fluctuations and 
solve the equations. We briefly conclude in section IV. 



II. PERTURBATION EQUATIONS 



We start with ten-dimensional dilaton gravity coupled to matter Lagrangian C m repre- 
senting brane winding and momentum modes: 

d 10 x e~ 2 * \R + 4(V0) 2 + e a4, C : 



S 



The field equations that follow from the above action can be found as 



iV + 2V M V,0 - -g^ [R + 4V 2 - 4(V0) 5 
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W - (a - 2)T, A V>, 



where T„ v is the matter energy-momentum tensor: 



(1) 

(2) 
(3) 
(4) 

v ~g dg^ V v »~""> ' (5) 
Eqns. (2) and (3) are obtained by directly varying the action with respect to the metric and 
the dilaton. Although for matter equations one in principle needs Lagrangian C m explicitly, 
in our case (4) can be deduced from the consistency of (2) and (3). 

It is well known that one particular combination of the field equations gives a constraint 
on initial data. Let to denote the unit normal vector field to an initial value hypersurface. 
One can easily show that (2) contracted with n^n v 



T = 



1 



d 



iV + 2V^0 - - 9lu , [R + 4V 2 - 4(V0) 2 ] - e°* T, v ) nV 







(6) 



does not contain any second order time derivatives. Firstly, it is well from the initial value 
formulation of general relativity that n^n u (R^ — g^R/2) involves only first order time 
derivatives of the metric (see e.g. [33]). On the other hand, the dangerous terms with two 
derivatives acting on dilaton in (6) can be reexpressed straightforwardly as 

V 2 + n^V M V,0 = /^V M V„0, 

= ^-(V)KV^), (7) 

where = — n^n u and = h^'SI v is the spatial covariant derivative along constant 
time hypersurface. It is clear that (7) contains only first order time derivatives and thus (6) 
can be viewed as a constraint on initial data. 
Assuming a metric of the following form 

rfs 2 = _ e 2A dt 2 + e 2B dx i dx i + e 2C ^2 ^ (g) 

(i = l,..,m) the energy-momentum tensor for a gas of branes winding over the compact 
manifold S p can be found as (see e.g. [22]) 

rp _ rp -VtlB , rp -mB-(j)+l)C 

- 1 - oo ± w'~- ~r-i m e , 

T Vj = 0, (9) 
T - - -T p~ mB ft , 4- — p -mB-( P +i)C x , 

A ah — i » c u ab \ c u abi 

p 

where T w and T m are constants corresponding to winding and momentum modes. In this 
paper we take S p to be Ricci flat. From (9) an equation of state p^ = uj^p can be deduced 
where u>i = 0, u a = —1 for winding modes and u>i = 0, uo a = 1 jp for momentum modes. Note 
that the equation for momentum modes is identical to radiation confined in the compact 
space. One can check from (5) and (9) that 

C m = -2p, (10) 

which is precisely the Lagrangian for hydrodynamical matter (see e.g. section 10.2 of [34]). 
In general we have V M T Ml/ ^ due to the coupling of C m to dilaton. 
To focus on D-branes one should set 

a = l, (11) 
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which is required to have the correct g s = dependence in the low energy effective action 
(1). Recall that we also have 

m + p = 9, (12) 

where m and p are the dimensions of the observed and the compact dimensions, respectively. 
In [24], we have found a particular solution to (2)-(4) 

ds 2 = -dt 2 + (t) 8 /( m + 3 ) dx l dx l + e 2C ° rfSj, 

e 4> = e <Po ^2(m-3)/(m+3) ) ( 13 ) 

e ^o = 4(p + 2)p £(p+1)Co = (p + 2)T m 

T„,(p+l)(m + 3) 2 ' pT w ' 

which has some appealing properties. For instance, it is possible to argue analytically 1 

and support the argument numerically that for arbitrary initial conditions (within a given 

ansatz) the fields approach (13) asymptotically. Moreover, for m = 3, the dilaton becomes a 

constant and thus the canonical Einstein and string frames become identical up to a constant 

scaling. 

Let us now consider linearized inhomogeneities on (13). We are interested in the case 
where the fluctuations are independent of compact coordinates 2 and choose the generalized 
longitudinal gauge (see [18]) in which the scalar metric perturbations only appear in diagonal 
entries. Therefore, the perturbed solution has the form: 

ds 2 = -e 2A {l + 25A)dt 2 + e 2B {l + 25B) dx l dx l + e 2C (l + 2dC) dT 2 , 

e* = e^(l + ^), (14) 

where A, B, C and <fi denote the background values depending only on t and 

6A = SA(t,x i ), 5B = 5B{t,x l ), dC = dC{t,x i ), &f) = &f){t,x l ). (15) 

Let us remind that in the longitudinal gauge 5A, 5B and dC coincide with diffeomorphism 
invariant quantities. 

In obtaining the linearized field equations we find it convenient to use (6) together with 
a particular combination of (2) and (3) 

+ 2V M V,0 = \g^p + e 4, T^. (16) 



1 As discussed in [30] and [31], it is possible to give analogous arguments for intersecting brane configurations 
with or without dilaton. 

2 This is indeed a good approximation since the extra dimensions are assumed to be very small. 
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As we pointed out above, (6) can be interpreted as a constraint on initial data . 3 Substituting 
(9) and (14) in these equations and keeping only the linear terms one gets 



25/> - m 50 -pdC + e~ 2B £ 8A kk + 2 [mB + mBB — 20] 54 + 

k 

mB - 20] 54 - 2mB5B = \pe' i> 5<p + \e*5p , (17) 
(m - 1)50; - 2h\) i - [(m - 1)0 - 20] 54; + pSCi - pBdd + 255/); = -e* +B 5T 0i , (18) 
50 + e- 2B [2(50;; - 54;; - p£7;; - (m - 2)50;;] - e' 2B ]T 50 fcfc + 2 



mi? — 



50- 



054 + p05i7 - 205> + 2 [200 - - m00] 54 = \pe <t> 5<p + ^e^5p , 
[25/) - 54 - pdC - (m - 2)50] . . = 0, i^j 



SC - e~ 2B 



]T ZCkk + [mB - 20] £7 = ^(50 + ^5p + e*T aa 5/) + e^5T aa , 



m[m 



- 1)0 - 2m0] 50 + [mp0 - 2p0] 80 + [4m00 - 400 - m(m - 1)00 



(19) 
(20) 
(21) 



54 + 



40 - 2m0] 5> + e- 2B ]T [2<50 fefe - (m - l)50 fefc - p80 kk ] = pe^&j) + e^p , 



(22) 



where (17)-(21) follow from (16), and (22) can be obtained from (6). In the above equa- 
tions the dot and the sub index i denote partial differentiations with respect to t and x\ 
respectively, and we have not used the summation convention. 

Note that in (18) we introduce a non-diagonal component for the linearized energy mo- 
mentum tensor ST 0i . Although for long wavelength fluctuations this equation is satisfied 
identically, it is in general necessary to add 5To; since the perturbed metric functions de- 
pend on both t and x\ which gives non-zero Ra (indeed (18) is the (Oi) component of (16)). 
Moreover, one can see that without ST 0i the matter equations constrain the metric functions 
too much, see (24) below. The appearance of this non-diagonal term for winding modes can 
be deduced from the Born-Infeld (BI) action as follows. When the metric functions depend 
only on t, a static brane winding S p and located at a constant position becomes 
an extremum of the BI action (i.e. the embedding coordinates become harmonic maps). 
Using this solution in the energy-momentum tensor obtained from BI action and smearing 
the delta function singularity by taking a continuum average for a gas of branes, one gets 
the winding part of (9) [9] . However, when the metric is perturbed with t and x % dependent 
functions, the static brane solution should be modified x % = x l + 5x % (t) and this correction 



It is possible to show that (6) and (16) are equivalent to (2) and (3) with a = 1. 
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induces non-zero ST 0i . 

Linearizing the matter equation (4) we find that v = a component is identically satisfied 
while v — and v — i equations give 



Sp + (mB + P C)5p + (mSB + P SC)p + SCT a a + C ST a a + e A ~ B STj = 0, (23) 



Of 



m + l)B + pC - (j) ST 0l + e A ~ B [ddT a a + pSPi - p5A t ] = 0. (24) 



Eq. (23) has actually two separate parts corresponding to winding and momentum modes. 
In obtaining Sp and ST^ from (9), one should also consider the variations ST W and ST m . 
Therefore, there are seven variables which are ST W , ST m , STq{, SA, SB, SC and Sp. Equations 
(23) and (24) determine the evolution of source fields ST W , ST m and ST 0i , which in turn 
should be substituted in (17)- (22). It is possible to solve (20) as 

2&f>-5A-pdC-(m-2)SB = 0, (25) 

which reduces the number of unknown functions by one. On the other hand, (18) and (22) 
are not dynamical evolution equations since they only contain first order time derivatives 
and thus merely restrict the solution space. Ignoring them, it is easy to see that the number 
of equations and unknowns are the same. 



III. LONG WAVELENGTH FLUCTUATIONS 



Let us now focus on long wavelength perturbations. As pointed out in [18] an instability 
on such scales would be a problem for the stabilization mechanism. In this regime all spatial 
derivatives can be ignored in field equations which in turn implies that (18) and (24) are 
identically satisfied. On the other hand, from (23) one obtains 

Sf w = Sf m = 0. (26) 

Since constant shifts of T w or T m alter the solution (13), they are not honest perturbations 
and one should set ST W = ST m = 0. In other words, (26) gives the zero modes which must 
be ignored. Solving Sp from (25) 4 and using the background values from (13), we find that 

4 For long wavelength fluctuations (25) can be viewed as the residual gauge symmetry fixing. 
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(17), (19), (21) and (22) respectively become 

m - 2m 4- 12 ^ _ 8m5i3 _ 1q p 6A + 2 ^ m + = n r 2 ?^ 

(m + 3)f (m + 3)t (m + 3) 2 t 2 (m + 3) 2 t 2 ' 1 J 
ar_ 8 ^ 20(50 10p&4 2p(m + 2) (SB _ 

(m + 3)t (m + 3)* (m + 3) 2 t 2 (m + 3) 2 t 2 ' 1 j 

£ + + ^ + ^ = , (29) 

(m + 3)i (m + 3) 2 t 2 ' v ; 

12&4 4(m - 6)5B 20p5A Ap(m + 2)50 



(m + 3)f (m + 3)f (m + 3) 2 t 2 (m + 3) 2 t 2 
Note that 5C decouples from 8A and 5B. As a consistency check of these equations, we find 
that the time derivative of the constraint (30) is proportional to four times (27) plus m times 
(28) and thus identically satisfied, as it should be. This does not mean however that (30) 
is empty. Rather, one should first solve (27)- (28) and then use (30) to narrow the solution 
space. 

Eq. (29) can easily be solved for dC which gives 

c^y/,\ cos(alnt) sin(alnt) 

SC(t) = ci K fb + c 2 K tb \ (31) 



where a, b are positive numbers depending on m (e.g. for m — 3, a — y 29/12, b — 1/2). On 
the other hand, the coupled equations (27) and (28) give 

5A= (m + 2) h + k 2 t-( 7 + m )/(«+3) + ^ ( m _ 2 ) t (-9+m)/(m+3) + ^ (39 + Um ) f 2 

(SB = 5fci - k 2 t-(7+™)/(™+3) + ^ t (-9+m)/(m+3) + ^^3 + m ) j2 (33) 

where fci,..,&4 are arbitrary constants. Substituting these solutions into (30), we observe 
that it yields 

h = 0, (33) 

and thus the growing mode in (32) is eliminated by the constraint equation. As a result, we 
find that there are no instabilities to this order. 



IV. CONCLUSIONS 



In this paper, we consider long wavelength fluctuations on the background of a D-brane 
gas in ten dimensional dilaton gravity constructed in [24] . The most attractive feature of the 
solution is that it yields stabilized internal dimensions. Moreover, when the observed space 
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is three dimensional the dilaton also becomes a constant. The main result of this paper is 
that the long wavelength perturbation of the internal space performs damped oscillations. 
Other fluctuations, including the dilaton, are found to decay to a constant. Therefore, the 
stabilization mechanism based on D-brane winding and momentum modes remains valid in 
the presence of linearized inhomogeneities. 

Technically, the perturbation analysis for D-branes is similar to the string gases studied 
in [18]. The main difference is that in our case there appears to be a growing mode which is 
eliminated by the constraint equation. Although the presence of this mode would not affect 
the internal space, it would ruin the stabilization of dilaton for m — 3, which is crucial in 
the effective field theory obtained after compactification. 

Higher dimensional D-brane gases offer an analogous scenario to string gases for the 
stabilization of the extra dimensions. Moreover, D6 branes happen to fix dilaton in three 
dimensions. This mechanism can be useful in topologically non-trivial compactifications 
where the winding strings do not appear in the spectrum or they are not capable of fixing 
volume modulus. Of course, the toy model we consider is immature and it should be devel- 
oped in many different directions. However, the results in BGC show that string theory has 
the right ingredients to solve the stabilization problem of extra dimensions in cosmology. 
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